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ASYMPTOTICS OF SOLUTIONS TO SEMILINEAR 
STOCHASTIC WAVE EQUATIONS 

By Pao-Liu Chow^ 

Wayne State University 

Large-time asymptotic properties of solutions to a class of semi- 
linear stochastic wave equations with damping in a bounded domain 
are considered. First an energy inequality and the exponential bound 
for a linear stochastic equation are established. Under appropriate 
conditions, the existence theorem for a unique global solution is given. 
Next the questions of bounded solutions and the exponential stability 
of an equilibrium solution, in mean-square and the almost sure sense, 
are studied. Then, under some sufficient conditions, the existence of 
a unique invariant measure is proved. Two examples are presented to 
illustrate some applications of the theorems. 

1. Introduction. Semilinear stochastic wave equations arise as mathe- 
matical models to describe nonlinear vibration or wave propagation in a 
randomly excited continuous medium. To be specific, the equation may take 
the form 



(1.1) d'^u{x, t) = (?I^u - 2a dtu{x, t) + f{u) + a{u)W{x, t) 

in a bounded domain T) in R'^, subject to some homogeneous boundary 
and initial conditions to be specified later. Here 9i = ^, A is the Laplacian 
operator, and c and 2a are some positive constants known as the wave speed 
and the damping coefficient, respectively. The nonlinear functions / and 
a are given, and W{x^t) = dtW{x,t) is a spatially dependent white noise, 
where W{x,t) is a Wiener random field. In previous papers [3, 4], we studied 
the local and global solutions of this type of equation without damping [a = 
0), where the nonlinear terms / and a may admit a polynomial growth. As 
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a sequel to our previous work, this paper is concerned with some quahtative 
asymptotic behavior of solutions to equation (1.1) in a bounded domain D 
as t ^ oo. In addition to the global existence of solutions, we are interested 
in the questions of boundedness, asymptotic stability and the existence of 
a stationary solution or an invariant measure. For a solution of the wave 
equation to reach a statistical equilibrium, it is imperative to include the 
damping term in equation (1.1) so that, in the physical term, the fluctuation- 
dissipation principle may hold. As a simple example, consider the randomly 
perturbed wave equation in one dimension, 

d^u = c'^dlu-2adtu + W{x,t), t>0, x G 2? = (0, vr), 

(1.2) 

u{x,0) = h{x), dtu{x, 0) = 0; u{0,t) = u{7r,t) = 0, 

where /i is a given continuous function and the Wiener field W is assumed 
to have the Fourier series representation 

oo 
n=\ 

where {6n(i)} is a sequence of independent copies of standard Brownian 
motions in one dimension, {(7„} is a sequence of reals such that Y^=\ < 
oo and (/>n = ^2/7rsin nx, n= 1,2,..., are the normalized eigenfunctions 
associated with the problem (1.2). Then, by means of the eigenfunction 
expansion, (1.2) can be formally solved in the case c > a to give 

oo 

(1.3) u{x,t) = ^Un{t)(t)n{x), 

n=l 

where 

(1.4) M„(t) = /i„e~"*cosu;„,t + — / e""^*"*) sinu;„(t - s) d6„(s) 

with hn = h{x)(j)n{x) dx and ujn = \/ {ncy — for n = 1, 2, . . . . By some 
simple calculations, we obtain the mean 

EUn{t) = hne~"^ COSLOnt 

and the variance 

VarK(t)} = /'*e-2-sin2(a;„s)ds^ 

\uJnJ Jo 4a 

as t ^ oo. Then it follows from (1.4) that the solution u{x,t) is a Gaussian 
random field with mean Eu{x,t) 0, and covariance function 

lim Cov {u{x,t),u{y,t)} = 7— — (t>n{x)cl)n{y)- 
t-*oo 4:a\nc/ 
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In fact it can be shown that the solution u{-,t) converges in the mean- 
square to a Gaussian random field u(-) with the above covariance function 
and its probability law is the invariant measure for equation (1.2). On the 
other hand, without the damping (a = 0), we would have Eun(t) = hn cos net 
and 

Var{n„(t)} = {l/2){an/ncf[t - (l/2?ic) sin2nct] ^ oo as t ^ oo. 

So the asymptotic solution will cease to exist. Clearly this can also happen in 
the nonlinear case. In fact it was shown that, with a cubic nonlinearity, the 
solution may explode in finite time [3] unless there exists a certain energy 
bound. As to be seen, the dissipation and the energy bound for a semilinear 
wave equation such as (1.1) are two major ingredients to ensure proper 
asymptotic behavior of its solutions. 

Our initial work on semilinear stochastic wave equations [3] was stim- 
ulated by two interesting papers by Mueller [18, 19] on the existence of 
large-time solutions to some nonlinear heat and wave equations with noise. 
If such a solution exists, it is natural to investigate its asymptotic behav- 
ior as t ^ oo. By a semigroup approach, asymptotic solutions to semilin- 
ear stochastic evolution equations have been studied by many authors. For 
the problems of boundedness and stability, see, for example, the papers 
[2, 10, 12, 13], and for the existence of invariant measures, we mention the 
articles [5, 8, 14, 16], and the book [7] for further references. In concrete 
terms, most of the above-mentioned results are applicable to the parabolic 
or dissipative type of stochastic partial differential equations. The asymp- 
totic solution of a stochastic hyperbolic or wavelike equation was studied in 
[15] by the method of averaging. To our knowledge the asymptotic solutions 
of the semilinear wave equations under consideration have not been treated 
in the literature. For the deterministic case, the analysis of hyperbolic equa- 
tions relies heavily on the so-called energy method ([23], [25], Chapter 4). 
Therefore the associated energy function plays an important role in the 
asymptotic analysis. Similarly we shall adopt the stochastic version of the 
energy method in the current study. In fact, to obtain the crucial exponential 
estimates, it is necessary to introduce a pseudo energy function, which can 
be interpreted physically as adding an artificial damping to the system. For 
some related works on stochastic wave equations, we mention the interesting 
papers [17, 20, 21], among many others. 

2. Summary of results. In Section 3 we present three technical lemmas. 
In Lemma 3.1 we prove the existence of a unique solution and the energy 
equation for a linear stochastic wave equation. By introducing a pseudo 
energy function, a key exponential estimate is established in Lemma 3.2. 
Then it is shown in Lemma 3.3 that the pseudo energy function is equivalent 
to the usual energy function. 



4 



P.-L. CHOW 



The global solution to a class of semilinear stochastic wave equations 
of the form (4.1) is treated in Section 4. Under the locally bounded, local 
Lipschitz conditions in the Sobolev space and an energy inequality given 
by conditions (Al)-(A4), the results of the existence and uniqueness of a 
global solution are stated and proved in Theorem 4.1. The proof is based on 
a smooth if ^-truncation technique and some probabilistic inequalities. 

In Section 5 we consider the boundedness of solutions in a mean-square 
sense as t ^ oo. Assuming that, in addition to Conditions A, the nonlinear 
terms satisfy a set of growth conditions (B1)-(B3), Theorem 5.1 shows that 
the solution is bounded in mean-square, while with a slightly stronger as- 
sumption, it is proved in Theorem 5.2 that the solution is ultimately bounded 
in the mean-square sense. 

Then some questions of the asymptotic stability of the null solution are 
considered in Section 6. Under Conditions B with an exponential integrabil- 
ity condition on the parametric functions 9{t) and p{t), Theorem 6.1 shows 
that the null solution is asymptotically, exponentially stable in mean-square. 
If = p = 0, as stated and proved in Theorem 6.3, the null solution becomes 
exponentially stable almost surely. 

So far the stochastic wave equations under consideration admit that non- 
linear terms satisfy only a local Lipschitz condition in the space H^. In 
particular, the nonlinear terms are allowed to have a polynomial growth. 
For the existence of an invariant measure, this poses a challenging open 
problem as yet to be resolved. Even in the case of globally Lipschitzian non- 
linearity, the existence result does not seem to have been proven. Hence, in 
Section 7, we shall prove an existence theorem (Theorem 7.1) by assuming 
that the nonlinear terms are globally Lipschitzian and have linear growth. 
Technically our proof follows the approach of Da Prato and Zabczyk [7] 
by adapting their method for the strongly dissipative equation, such as a 
parabolic equation, to our hyperbolic problem. Finally two examples are 
provided in Section 8 to illustrate some applications of our theorems. 

3. Energy equation and exponential estimate. Let V C R'^ be a bounded 
domain with a smooth, say, boundary dT>. We set H := LF'iT') with 
the inner product and norm denoted by (•,•) and || • ||, respectively. Let 
= W^''^{'D) be the Sobolev space of order k with norm || • and 
denote by Hq the closure in of the set of all functions with compact 
support in V. The dual space of is given by [1]. 

Let P) be a complete probability space for which a filtration J^t of 

sub-c7-fields of T is given. Let W{x, t), x G P, t > 0, be a continuous Wiener 
random field defined in this space with W{x,0) = 0. It has a zero mean, 
EW{x,t) = and covariance 



(3.1) E[W{x,t)W{y,s)] = itAs)r{x,y), x,yeV, 
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where (t A s) = min(i, s) for <t,s <T and the covariance function r{x, y) 
is bounded so that 

(3.2) snpr{x,x) < vq. 

Let (t(x, t) = a{x, t, to) for t > 0, x G 2? and a; G be a continuous .7"t-pre- 
dictable random field satisfying the condition 

(3.3) E r \\a{-,t)\\Pdt<oo 

Jo 

for p>2. Then it can be shown that the stochastic integral 

(3.4) M{x,t)= [ a{x,s)W{x,ds), t>0, xeV, 

Jo 

is well defined and = M{-,t) is a continuous ff-valued J-t martingale 
(see the Appendix). It has mean EM{x,t) = and covariation operator Qt 
defined by 

{{M.,g),{M.,h))^= [\Qsg,h)ds 
Jo 

for any g,h& H, where the kernel function q{x,y,t) of Qt, defined by 

{Qt9){x)= I q{x,y,t)g{y)dy, 
Jv 

is given by 

Q{x,y,t) =r{x,y)a{x,t)a{y,t). 

In view of conditions (3.2) and (3.3), it can be shown that (see the proof of 
Theorem A.l) 

E\\Mtr<CpiT)E r \\ai;tWdt 
Jo 

for some positive constant Cp{T). It is worth noting that the stochastic 
integration in (3.4) is taken with respect to an L^-bounded integrand at, 
instead of a Hilbert-Schmidt operator-valued process as usually done (see 
[6], Chapter 4). This version of stochastic integral will be needed later on to 
deal with equations with pointwise (in x) multiplicative noises (see Example 
1). Since we have not been able to find a reference for this type of integral, 
it will be defined in the Appendix. 

Now we consider the initial boundary value problem for the linear damped 
hyperbolic equation with a random perturbation, 

[dt + 2adt - A{x, D)]u{x, t) = f{x, t) + dtM{x, t), 0<t<T, 

(3.5) u{x,0) = uq{x), dtu{x,0) = vo{x), xGV, 

u{-,t)\dv = 0, 
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where a is a positive parameter, D = dx denotes the gradient operator and 
A{x,D) is a strongly ehiptic operator of second order of the form 

d 

(3.6) A{x,D)^{x) = J2 d,^[a'\x)dx^^ix)]-b{x)^ix). 

In addition, the coefficients a*-' = a-'* and b are assumed to be smooth func- 
tions that satisfy 

d 

for some constants ai > qq > 0. 

To consider (3.5) as an Ito equation in a Hilbert space, we set ut = 
u{-,t),vt =v{-,t) and so on, and rewrite it as 

dut = vtdt, 

(3.7) dvt = [Aut - 2avt + ft] dt + dMt , < t < T, 

^^0 = 9, vo = h, 

where the domain 'D{A) = n Hq, g G H^,h G H and Mt is regarded as 
an i?-valued Wiener martingale. Condition (3.6) implies that {—A) is a 
self-adjoint, strictly positive linear operator in H = L'^{D) and its square 
root B = \J —A is also a self-adjoint, strictly positive operator with domain 
T>{B), which is a Hilbert space under the inner product {g,h)B '■= iBg,Bh) 
and norm ||5'||_b = 

[Bg^BgY'"^ (see [24], Chapter 1). Since D{B) = Hi, for 
convenience, we define || • ||i = || • ||b in the subsequent analysis. As usual, the 
Ito differential equation (3.7) is interpreted as a stochastic integral equation: 

ut = uo+ / Vsds, 

vt = vo+ / Ausds — 2a / Vsds+ fsds + Mt. 
Jo Jo Jo 

Introduce the Hilbert space H = {H^ x H) with Ho = {Hq x H), equipped 
with the norm defined by 

\mn = {\\u\\l + \\v\\Y^' = {\\Bnf + \\vf}'/' 

for any (f> = {u;v) G H. Let TL* = {H~^ x H) denote the dual space of H. 
Define the energy function e(-) : TC — > R"*" = [0, oo) as 

(3.9) e{4>):=e{u;v) = \\Buf+\\vf for (f) = {u; v) e x H . 



Notice that the norm Hf/)]!?^ = \/^{4') is ^-Iso called an energy norm. In what 
follows, we denote the "H-norm || • simply by || • || when there is no con- 
fusion. 



SEMILINEAR STOCHASTIC WAVE EQUATIONS 



7 



Now by regarding (3.8) as a stochastic evolution equation in 7i* in the 
distributional sense, we have the following lemma: 



Lemma 3.1 (Energy equation). For <^o = {^o^vq) G 7i, let ft he a con- 
tinuous predictable process in H and let Mf be a continuous H -valued mar- 
tingale with covariation operator Qt such that 

(3.10) eI^J^ Wftf dt + j\rQtdt'^ < oo, 

where Tr denotes the trace operator in H. Then equation (3.8) or (3.7) has a 
unique solution (j)t = {ut]Vt) which is a continuous TC-valued semimartingale. 
Moreover, it satisfies the energy equation 

e{4>t) =e{(f)o) -4a f \\vs\\'^ds + 2 [ {vs,fs)ds 

+ 2 / {vs,dMs)+ [ Tr Qsds a.s. 
Jo Jo 

for t G [0,r], where the energy function e(-) on TC is defined by (3.9). More- 
over, the inequality 

(3.12) E sup e{(l)t) < Ci + C2E [ + TYQJds 

t<T Jo 

holds, where the constants Ci,C2 depend on p,T and the initial conditions. 



Proof. Since the idea of the proof is similar to that of Lemma 2.1 in [4] 
with the Laplacian replaced by A, we will only sketch the proof. The only 
difference is that, instead of using Friedrichs' mollifying approximation, we 
adopt a finite-dimensional projection. 

To this end, since A is strongly elliptic and self-adjoint, it has a complete 
orthonormal set of eigenfunctions {<^n} with corresponding eigenvalues {A^}. 
Let Pn'.H ^ Hn be defined by Pnh = X)fc=i(^) ^k)v'k, where Hn is a finite- 
dimensional subspace of spanned by {ipi, . . . ,ipn}- Apply the projector 
Pn to equation (3.7) to get 

du'l = v1 dt, 

(3.13) dvl' = [Au'l -lav"^ + fJ'jdt + dM]^, 0<t<T, 

„.n _ n n _ <n 

— g , Vq — n , 

where we set = PnUt, . . . , /i" = Pnh. The finite-dimensional linear system 
has a unique ^f-adapted continuous solution (ti";u") in {Hn x Hn) C {Hq x 
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H). In particular, by the Ito formula, the following energy equation holds: 



(3.14) 



e(nr,<) = e(g",/i")-4a t\\v^fds + 2 f\v:,f^)ds 

Jo Jo 



+ 2 f\v^,dM2)+ TTrQ^ds 
Jo Jo 



a.s. 



By means of the simple inequality 2(f", /") < (a||f"|p + — H/'^lP) and the 
B-D-G (Burkholder-Davis-Gundy) inequality ([6], page 82), 



ft 

E sup 

0<t<T 



rt ( fT ) 1/2 

! J «,dM;) <8^|y \K\\ TtQ'^dsj 

<8e{ sup IK II /^IKIlTVQ^ds ^ 1 

lo<t<T Jo ) 

rT 

<\E sup ||urf + 36/ TrQ^ds, 

0<i<T Jo 

we can deduce from (3.14) that 

(3.15) E sup e(ur,0<Cii? All/.f + Trg:}(is 

0<t<T Jo 

for some constant Ci > 0. Let X = L'^{Q.;C{[Q\T],HI x i?)) with norm 
||(ii;?;)||x = {£^supQ<j<y e(nf,ft)}^/^. Then X is known to be a separa- 
ble, reflexive Banach space ([22], page 218). In view of (3.15), the sequence 
{{u^;v'^)} is bounded in X so that there exists a subsequence {(«"''=; f"*=)} 
that converges weakly to {u]v) E X. 

In fact, we can show that the subsequence converges strongly in X. To 
do so, denote the subsequence again by {(n";f")} and set (^»"";u™") = 
(n™;t;™) — (n";!;"). It suffices to show that {(u";f")} is a Cauchy sequence 
in X so that || («'"'"; u™")||x ^ as m,n^oo. In view of (3.13) and (3.14), 
the difference sequence satisfies the energy equation 

e(uP;t;P)=e((7"''^;/i"''^)-4a C Wv^J'^'f ds + 2 f\vT'' , fr"") ds 

(^•16) . Jo ^ Jo 

+ 2 f {vf'',dM^'')+ f TVQ^^ds a.s. 
Jo Jo 

By similar estimates that lead to (3.15), we can obtain 

E sup e(nP,?;r")<C2(e(5"''^,/i'"")+^ A||/f"||2 + TrQ™"}(i4 

0<t<T I Jo J 

(3.17) 
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for some constant C2 > 0. Since the right-hand side of (3.17) tends to zero 
as m,n^oo, it follows that \\{u"^;v"^) — {u"';v"')\\ -^0. Hence {(n";ti"')} is 
a Cauchy sequence in X and hm„^oo('«"; 'w"') = iu',v) strongly as claimed. 
Due to this strong convergence, it is easy to show that the limit (u; v) is the 
unique strong solution with the depicted regularity. Moreover, we can take 
the limits termwise in (3.14) to obtain the energy equation (3.11). Then the 
energy inequality (3.12) follows easily. □ 



Notice that, due to the lack of required smoothness of solutions, the gen- 
eral Ito formula does not hold here. As in the deterministic case, the energy 
equation and the associated inequalities are the key to proving the existence 
and regularity results for stochastic hyperbolic equations. 

Owing to the dissipation term in (3.11), in contrast to the energy in- 
equality (3.12), it is possible to obtain an exponential estimate for the mean 
energy. To this end, we introduce a pseudo energy function 

(3.18) e^{(f)):=e^{u;v) = \\Buf + \\v + Xuf lor ueH^,veH, 
where A > is a parameter. Let = v + An. Then we can write 

(3.19) e^{u;v) = e{u;v^) = e{u;v) + 2X{u,v) + X^\\uf . 

Since A is strongly elliptic and strictly positive, its smallest eigenvalue 771 
can be characterized as ([25], page 62) 

(3.20) 7?i= inf ^^>0. 

geH\gj^o \\g\\^ 

Lemma 3.2 (Exponential estimate). Let the conditions for Lemma 3.1 
be satisfied such that (3.10) holds for any T > 0. Then if 

(3.21) A<Ao:=min|- 



, 2 ' 4a , 

there exists ai G (0, A) such that the following inequality holds: 

(3.22) i?e\,^t) < e\<^o)e-"^* + f e-^^^^'"'^ E\ — \\f,f + T,q\ ds. 

Jo [ai J 

Proof. It follows from (3.7) that {ut\v^) satisfies the perturbed system 

dut = [vi — Xut] dt, 

(3.23) dv^ = [Aut + A(2a - X)ut - {2a - X)v^ + ft] dt + dMt, 

uo = 9, vo = h, <t <T. 
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By applying Lemma 3.1 to the above system and noting e'^(ii( ; w^) = e{ut;Vf), 
the pseudo energy function (3.18) satisfies 

de\ut;vt) = 2[X{2a - X){ut,v^) 

(3.24) - X\\ut\\l - {2a - X)\\v^f + ift,v^) + iTrQ^] dt 

+ 2iv^, dMt), 

with e'*'(uo;fo) =^{uq;vq). Now, in view of (3.20), we have, by using some 
simple inequalities. 



X{2a - X){u,v^) - X\\u\\l - (2a - X)\\v 



A||2 



< X{2a - X)^^\\v^\\ - X\\u\\l - (2a - A)||^;^ 



< X{2a - X) 



X Ikll? 1 II Aii2 



VI 4A' 



X\\u\\l - (2a- A)||?; 



A||2 



2aA\,, ,,9 3a,, \,,r, A,, ,,9 3A,, \,,9 
< -A( 1 ]\\u\\l - —\\v^f < --\\u\\l - — . 



r/i / ■■ - 4 " " - 2 
The above result together with the fact that 



{fuv')<^\\vY + \\\ftf 



imply that 



(3.25) 



A(2a - A)(n, v^) - X\\u\\i - (2a - X)\yr + ift,v^) 

<-^{Ml + \\v' 



In view of (3.25), equation (3.24) yields 



(3.26) de^{ut;vt)<-Xe^{ut;vt)dt + 



dt + 2{v^,dMt), 



which can be integrated to get the desired inequality (3.22), after taking the 
expectation, with any ai < A. □ 

It is easy to show that the energy norms induced by e and are equiv- 
alent. In fact, the following lemma holds. 



Lemma 3.3. For any X G (0,/ii), the inequality 



(3.27) 



m + x 



Km - X 



e[u; V) 
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holds, where fii = (■\/4r/i + A^). Moreover, we have 

(3.28) Ee{(t>t) < A'(A)|e((Ao)e-"i* + ^{^Jfsf + T>Q.) ds^ 

where K{\) = {i^ii + A)/ (/ii — A) . 

Proof. By definition (3.19), 

e^{u;v) =e{u;v) +2X{u,v) + A^||ti||^. 
It follows that, for any /3 > 0, 

e\u; v) < e(n; v) + {l + 0) j \\u\\l + -\\vf 

by choosing (3 = \{\/ {Arji/X'^) + 1 - 1}. 
On the other hand, for any 7 > 0, 

e^{u;v)>e{u;v)- |(7 " 1) Ikll? + ^Ikf | 



by taking -y = ^{Vl + 4r]i/X^ + 1}. 

Therefore, we have verified (3.27), and the result (3.28) is now an direct 
consequence of (3.22) and (3.27). □ 

4. Semilinear stochastic hyperbolic equations. Let us consider the initial 
boundary value problem for the hyperbolic equation 

d^u{x,t) = [A{x,D) - 2adt]u{x,t) 

+ f{u, Du, X, t) + cj(n, Du, x, t) dtW{x, t), t> 0, 

(4.1) 

u{x,0) = uo{x), dtu{x,0) =vo{x), xeVcIV^, 
u{-,t)\dv = 0, 

where, in contrast to the linear problem (3.1), f{s,y,x,t) and a{s,y,x,t) for 
xSP, t>0, sSR and y G R"^ are continuous functions, and Wt = W{-,t) 
is a continuous Wiener random field with covariance operator R with kernel 
r(x, y) for x,y gV. 

Similarly we rewrite the linear case as a system of ltd equations in 7i*: 

dut = vtdt, 

(4.2) 

dvt = [Aut - 2avt + Ft{ut)] dt + dMt{u) 
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or 

Ut = no + / Us ds, 

(4.3) 



Vt = Vo- 



I [Aus-2aVs + Fs{us)]ds + Mt{u), 
Jo 



where we set Ft{u) := f{u, Du,-,t), 
(4.4) Mt{u)= I 



and 5]t(-) : H is defined by T,t{u){x) := a[u{x), Du{x),x, t] for any u £ 

H^,xeV. 

We are interested in the large-time solutions of (4.1) when the nonlinear 
terms allow polynomial growth and are locally Lipschitz continuous. For the 
existence of solutions, we shall impose a set of sufficient conditions. In what 
follows, for r, s G R, let b{r) and k(r, s) be real-valued functions which are 
positive, locally bounded and monotonically increasing in each variable. Let 
us introduce a positive function 0(-;-):'H^R+ which is continuous and 
locally bounded, and let 

(4.5) e{u;v) < Q{u;v) < e{u;v) + c\\u\\i 

for any (n; v) £ 7i and some constants c > and /c > 2. As a shorthand 
notation, we set 

||St(n)||| = Trgt(n)= / r(x,x)[St(n)(x)]2dx 

Jv 

and impose the following conditions, which will be referred to later as Con- 
ditions A: 

(Al) A : n Hq H \s an elliptic operator as given in (3.6). 
(A2) Ft{-) -.H^ and St(-) :H'^ ^ H are continuous in t > 0. There exist 
functions b{r) and k{r,s) as indicated above such that, for any t > 

\\Ft{u)f + l\\j:tiu)\\l<bi\\uh) + q{t) 

for some locally bounded function q £ L^(R~^). 
(A3) In addition, 

\\Ft{n) - Ft{u')f + i||St(n) - M^^')fn < k{\\u\\i, \\u'h)\\u - u'\\l 

for any u, u' £ Hi,t > 0. 
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(A4) There exists a positive function depicted as above and constants 
Cj > 0, i = 1, 2, 3, and k < 1 such that 

\2iFs{Us),Vs) + \\^s{Us)\\ji}ds 



<ci+C2 / @{us;vs)ds - C3@{ut;vt) + Ke{ut;vt) 
Jo 

for any u. eC{n+;H^) nC^{n+;H) with vt = dtUt. 

Theorem 4.1. Let Conditions A hold true. Then, for uq = g £ Hi and 
Vo = h£ H , the problem (4.1) or the system (4.2) has a unique continuous 
solution u. G C{[0,T]; Hi) with dtu. € C{[0,T]; H) for any T > 0. Moreover, 
the following energy equation holds: 

e{ut,vt) = e{uo,vo) + 2 [ [{vg, Fs{us)) - 2a\\vsf]ds 

(4.6) . ^° 

+ 2 ( {vs,^s{us)dWs)+ I ||S,(u,)|||ds a.s. 
Jo Jo 

Proof. The existence proof is similar to that of Theorems 4.1 and 4.2 
in [3], and will only be sketched in steps as follows: 

Step 1. Hi-Lipschitz truncation. For iV > 1, let r7Ar(-):R"'" = [0,oo) — > 
R"*" be a Cq° function such that 



(4.7) m{s) = [ 



1, forO<s<A^/2, 
0, for s>N, 



and <r/Ar(s) < 1 for N/2 <s<N. For {u;v) eU, define Snu = r]N{\\u\\i)u, 
Fl^ (u) = r]Ni\\u\\i)Ft{SNu) and (n) = r/Ar(||u||i)St(57vn). Instead of (4.2), 
consider the truncated system 

dut = vtdt, 

(4.8) dvt = Aut dt + Fl^{ut) dt + (n) dWt, 

uq = g, vo = h. 

Step 2. Local solutions. By conditions (A2) and (A3), it can be shown 
that 

(4.9) \\Ft^ {u)f =m{\\u\\i)Ft{SNu) < ai{N) 
and 

(4.10) ||F,^(n) - Ff{u')f < a^iNMu - n')||i 
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for any u, u' £ Hi, v,v' £ H and for some positive constants ai, 02 depending 
on A^. Similarly, we can deduce that 

(4.11) \\^^{n)\\l<as{N) 
and 

(4.12) ||Ef (n) - J:^{u')\\% < a,{N)\\J{u - u')f 

for any u,u' € Hi, v,v' G H, where 0:3,04 are some positive constants de- 
pending on A^. 

Therefore the truncated system (4.8) satisfies the usual linear growth and 
the global Lipschitz condition. By invoking a standard existence theorem 
([6], Theorem 7.4), equation (4.3) has a unique solution (j)^ = {u^;v^) G 
L\n;Ci[0,T];H^xH)). 

Introduce a stopping time tn defined by 

T^=inf{t>0:||nf 111 >iV/2}. 

Then, for t < tn, ut = is the solution of (4.1) with dtUt = . As tn is 
increasing in N, let Tqo = lim7v->oo tn- Define ut for i < Too A T by ut = 
a t <TN <T. Then ut is the unique local continuous solution. 

Step 3. Global solutions. Assume condition (A4) is also satisfied. By 
taking the expectation, the energy equation reads 



Ee{utATN;vtATN) =e{uo;vo) + 2E / {vs, Fs{us)) ds 

Jo 

+ 2E {vs,T.s{us)dWs)+E ||S,(u,)||2jds. 

JO JO 

Letting pAr(t) = Ee{utATjs,]VtATN) ^-i^d invoking condition (A4), the above 
yields 

(4.13) PN{t) <[e{uo;vo) +ci]+C2 J pn{s) ds + KpN{t). 

Since k < 1, there exists C3 > such that 

Pn(.T) < C3[e{uo;vo) + ciY^^ = Ct- 
On the other hand, we have 

Pn(T) = Eex{uTATN) > E{I{tn < T)ex{uTATN)} 

> CE{\\uTArjll{rN<T)} > cf^YplTN < T} , 



2 

where I is the indicator function and C > is a constant. The above in- 
equality gives 

P{tn <T}< 4pn{T)/CN^ < ACt/CN'^. 
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Since the series J2'n=i P{'''n ^ T} converges, by the Borel-Cantelh lemma, 
we can conclude that 

P{too <T} = 

or Too > T a.s. for any T > 0. Now we let = utf^T^^ and denote its limit 
limjv_,oo ""i^ still by ut- Then ut is the global solution as claimed. The energy 
equation (4.6) can be verified by taking the limits termwise, as — > oo, in 
the energy equation for the A^th truncated system (4.8). □ 

Remark. In the above theorem, for simplicity, we assumed that W{x, t) 
is a scalar Wiener random field. Under an obvious modification, Theorem 
4.1 and the subsequent theorems still hold true when W = (W^^\ . . . , W^^^) 
is a fc-vector- valued Wiener random field and a = (cJi, . . . ,0"^) is another k- 
vector- valued, predictable random field such that the product a{-)W{-) = 
T!j=i'^jW^^\-) is interpreted as a dot product. 

5. Bounded solutions. In view of (4.2), we rewrite the hyperbolic system 
(4.2) as 

dut =vtdt, 

(5.1) 

dvt = [Aut - 2avt + Ft{ut)\ dt + T.t{ut) dWt, t> 0, 

with a given initial state {uo;vo) which is an J^q random vector in 7i. For 
the existence of bounded solutions, we shall impose Conditions B as follows: 

(Bl) There exist $ G C^(i7^; R+) with Frechet derivative G C{H^;H) and 
p. G C(R+ X H^;H) such that Ft{u) = -i$'(u) +pt{u) for any u e 
and 

ci <^>(n)<C2(l + ||u||^) 

for some constants ci and C2 > 0, > 2. 
(B2) There exist constants (3i > and 7i,(5i G R with i = 1,2,3, and es- 
sentially bounded functions and p which are locally integrable such 
that 

{<^'{u),u) > (3i<^{u) - -fi\\u\\l - 6i, 
\\pt{u)f </32Hu)+-,2\\u\\l+e{t) 

and 

\\^t{u)fn<P3^u)+jMl+p{t) 
for any u G Hi and t > 0. 



16 P.-L. CHOW 

(B3) The above constants satisfy 

(/3i - \)\^ - /33A - /32 > 0, 

(71 - i)A^ + 73A + 72 <0. 

Theorem 5.1 (Bounded in mean-square). Suppose that Conditions A 
and B hold true. Given uq £ and vq & H being random variables such 
that 

E{e{uo;vo) + 'i>iuo)} < 00, 

then the solution of the problem (5.1) is bounded in mean-square. Moreover, 
there exist positive constants Ki and 02 > such that 

E{e{ut;vt) + <^{ut)} 

(5.2) < i^i|£;[e(no; ^;o) + ^(^o; ^^o)]e-"^* 

+ ^*e-"2(*-') je{s) + p{s) ds^ + 2\5i\ Vt>0. 

Proof. By applying Lemma 3.2 to (5.1), as in (3.24), we obtain the 
perturbed energy equation 

de\ut;vt) = 2[\{2a - \)iut,v^) - X\\ut\\l - {2a - \)\\v^f 

(5.3) +iFt{ut),v^) + l\\My't)fR]dt 
+ 2dM^{u;v), 

where we set dM^{u]v) = {v^,i:{ut)dWt). As in (3.27) in Lemma 3.2, for 
A < {a A r]i/Aa}, the above yields 

de\ut;vt)<-X{\\ut\\l + l\\v^f)dt 

(5.4) 

+ [2{Ft{ut),v^) + \\J:t{ut)\\l]dt + 2dMt^(u;v). 

By assumptions, 

which, in view of condition (B2), implies that 
(5.5) 
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Define a superenergy function J -.TC ^ TV^ by 
(5.6) J{u;v)=e{u;v) +^{u), 

with J-^ = e^ + <I>. By applying (5.4), (5.5) and condition (B2) to (5.3), we 
obtain 

dJ^{ut;vt) < -Xe^{ut;vt)dt-(^f3iX-^-l33^^ut)dt 

+ [^\\vtf + (71A + y + 73) WutWl + \e{t) + p{t) + A5i| 

+ 2dMi-^(ii;i;). 
By invoking condition (B3), the above inequality gives 

dJ^{ut;vt) < -^J^{ut;vt) dt + IjOit) + pit) + X5i\ dt 

(5.7) 

+ 2dMt^{u;v), 

which implies that 

Ej\ut;vt)<EJ^{uo;vo)e-^'/^ 

/(^) + Ks) 

for all t > 0. Since, by assumption, 6 and p € Lj'Qj,(R+) are essentially bounded, 
we have 

EJ^iut; vt) = E{e^{ut; vt) + $(nt; vt)} < 00. 

Now, by invoking Lemma 3.3, J{u; v) < CJ'^{u\ v) for some C > 0. Therefore, 
the result (5.2) holds with some constant Ki > and 0.2 = ^- □ 

In fact, under somewhat stronger assumptions, it is possible to show that 
the solution of (5.1) is ultimately bounded in mean-square, that is, 

limsup-E{||uf llf + ll^'tll^} < 00. 

Theorem 5.2. Assume that Conditions A and B hold true with 5i = 0, 
and p € L^(R''~). Then the solution (f)t = {ut]Vt) is ultimately bounded in 
mean-square such that 

E sup {e{ut;vt) + ^{ut)} 

0<t<T 

(5.9) 

< K2E{e{uo; vo) + $(no)} + K3 [ [e{s) + p{s)] ds 

Jo 

for some positive constants K2 and 1^3. 



+ 



* -A(t-s)/2 



ds + 2\5i \ < 00 
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Proof. In view of (5.7), it is clear that 



A 



J {ut;vt) + 2 /- {us;vs)ds 



<j\uo;vo)+ f 
Jo 



Hence 



A 

A /•* 



ds + 2M^{u]v). 



(5.10) 



and 



(5.11) 



EJ^{ut;vt) + - f EJ^{us;vs)ds 

^ JO 



<Ej\uo;vo)+ f 
Jo 



9is) + p{s) 



ds 



E sup J^{uuVt)<EJ^{uo]Vo)+ I 

0<t<T Jo . 



1 

A 



e{s) + p{s) 



ds 



+ 2E sup \Mt^{u;v)\. 

0<t<T 

By means of the B-D-G inequahty for a submartingale, we can deduce 
that 

rt 



E sup \M^{u;v)\ =E sup 

0<t<T 0<t<T 



{v^,^s{Us)dWs) 







(5.12) 



(pi ^ J 

<3^|y^ {R^s{us)v^,^s{us)v^)dsj 

<3e( sup \\v^\\\{ r \\i:sius)\\Rds 

lo<t<T ) UO 

<lE sup \\v^f + 9Er\\J:,iu,)\\lds. 

0<t<T Jo 



Results (5.8) and (5.12) and condition (B2) imply that 
E sup J^{ut;vt)<EJ^{uo;vo) + lE sup \\v^f 

0<t<T ^ 0<t<T 



t r 



+ 18E r W^sMWlds + [ 
Jo Jo 



<EJ^{uo;vo) + -E sup ||<f 
^ 0<f<T 

+ 18^ r[P3<^{us) + ^3\\us\\i]ds 
Jo 



is) + p{s) 



ds 
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+ 



ds. 



Therefore, there exist positive constants Ci,i = 1,2,3, such that 

rT 

E sup J^{ut\vt) < ciEJ^{uq\vq) + C2 I E sup J^{Ur'^Vr)ds 
0<t<T Jo 0<T<s 



T 



C3 / [e{s)+p{s)]ds. 





From this together with the bound (5.10) and the Gronwall lemma, we 
can infer that there exists a pair of positive constants k2 , such that 

E sup J^{ut-vt)<k2EJ^{uo;vo) + k3 [ [e{s) + p{s)]ds, 

0<t<T Jo 

which, by Lemma 3.3, leads to the desired inequality (5.9). □ 

6. Asymptotic stability of solutions. Suppose that the hyperbolic system 
(5.1) has an equilibrium solution u = u£ T>{A) with v = 0. By a translation 
via u, without loss of generahty, we may assume that {u; v) = (0; 0) is an 
equilibrium solution. We are interested in the asymptotic stability of the 
null solution in the following sense. 

Definitions. 

1. The null solution (p = {u;v) = (0;0) of (5.1) is said to be asymptotically 
stable in mean-square in if 3 5 > such that, for \\4>o\\ < 6, 

lim E\\4>tf = 0, 

t— >oo 

and it is exponentially stable in mean-square if there exist positive con- 
stants K(6) and z/ such that 

E\\(j)tf <K{6)e-''* Vt>0, 

where = \\u\\l + ||f |p. 

2. The null solution is said to be a.s. (almost surely) asymptotically stable 
if 

pjlim \\^t\\=o\ = l, 

and it is a.s. exponentially stable if there exist positive constants K2{5),i'2 
and a random time T{uj) > such that 

\\<pt\\ < K2{6)e-''^* yt>T, a.s. 
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Remark. In view of the above definitions, it is clear tfiat the exponential 
stability implies the asymptotic stability. 

To proceed we assume that -Fj(O) = and Sj(0) = for any t > so that 
(j) = {u; v) = (0; 0) is an equilibrium solution of equation (5.1). In the stability 
analysis [11], it is often assumed that the global solution exists in the first 
place. Hence we suppose, under suitable conditions such as Conditions A, 
that the equation has a unique global solution. 

Theorem 6.1 (Stability in mean-square). Suppose that Conditions B 
hold true with the following provisions: 

(1) <I>(0) =0 and ^>(ti) > ifu^O. 

(2) In condition (B2), 6i=0 and there exists ckq > such that 



oo 



dt = C <oo. 



Then the null solution of equation (5.1) is exponentially stable in mean- 
square. Moreover, if (pQ = {uo;vo) is an -measurable random variable in 
TL satisfying 

E{e{uQ;vQ) + $(uo; vq)} < oo, 

then the inequality 

(6.1) E{e{ut;vt) + Hut)} < E{e{uo;vo) + Huq) + Cje-^^* 
holds for any t>0, where 03 = (ao A 02)- 

The theorem follows immediately from (5.2) in Theorem 5.1 and the sim- 
ple fact that e-°2(*-^)0(s) ds < e"""* /~ e"o''0(s) ds. In fact it is possible 
to show that the null solution is a.s. exponentially stable. Before stating the 
next theorem, we need a lemma which is a simple consequence of Theorem 
5.2. 

Lemma 6.2. Under the conditions for Theorem 5.1, the solution (fit = 
{ut;vt) is ultimately bounded in mean-square such that 

E sup {e^{ut]vt) + ^{ut)} 

0<t<T 

(6.2) 



< KiE{e\uo; vq) + ^uo)} + K2 f [0{t) + p{t)] 

Jo 

for some constants Ki,K2 > 0. 



dt 



With the aid of Theorem 6.1 and Lemma 6.2, we can prove the following 
theorem. 
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Theorem 6.3 (Almost sure stability). Assume that all of the conditions 
for Theorem 5.1 hold true with 9 = p = 0. Then the null solution of (5.1) is 
exponentially stable almost surely. Moreover, there exist positive constants 
C, V and a random variable T{uj) > such that 

(6.3) e{uuVt) + ^{ut)<C{e{uo;vQ) + ^uo)}e-''* a.s. 
for any t> T. 

Proof. Owing to Lemma 3.3, instead of (6.3), it suffices to show that 

(6.4) {e^{ut;vt) + ^{ut)}<CoE{e\uo;vo) + <!>{uo)}e-'\ t > T, a.s. 

for some constant Cq > and for A satisfying (3.21). To this end, decompose 
R+ as R"*" = lJ^o[^5"' + l] ^'^d consider the solution {ut]Vt) for n < t < n + 1 . 
Following the steps leading to (5.8) in the proof of Theorem 5.1, it can be 
shown that 

(6.5) E sup J^{ut;vt)<EJ^{un;Vn) + 2E sup \M^{u)\, 

n<t<n+l 0<t<T 

where we recall that J'^{u;v) = e^{u;v) + <^{u). As in (5.12), we have 



E sup \M^{u■,v) 

n<t<n+l 



<3e\ sup \\v^\\]( r^\\^s{us)fRds] ^ 

in<t<n+l J iJn J 

<3(e sup \\v^f] ^ {e r^\\^s{us)\\Rds] 

I n<t<n+l ) 1 Jn ) 



<3{E sup "-^"2 



1/2 ( rn+1 >, 1/2 

(/33 V 73)/ Ej\us;vs)ds\ 

n<t<n+l J I Jn J 

By making use of (5.11), the above gives rise to the upper bound 
E sup \Mt^{u;v)\ 

n<t<n+l 

( /-n+l ^ 1/2 

(6.6) <3{Ej\uo;vo)}^CK{P3yi3)J e"''''^ ds 

<6{i^j\no;fo)}|^(/33V73)}'^'e-"^/'. 
By taking (6.1), (6.5) and (6.6) into account, we get 

(6.7) e\ sup J^{uuVt)\<CQ{Ej''{uQ-VQ)]e-^^'^ 

Kn<t<n+1 ) 
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for constant Cq > 0. 

Therefore, by using the Markov inequality and (6.7), 



P 



sup J\ut;vt) > CoEj\uo;vo)e-''^/^ 

n<t<n+l 



< 



CoE{J^{uo;vo)}e~^^/^ 

Since J2'n=i^~"'^^^ < follows from the Borel-Cantelli lemma that 

there exists a random number N{uj) > such that, for n > N, 



sup J^{ut;vt) <Co{EJ^{uo;vo)}e 

n<t<n+l 



-nX/8 



a.s., 



which, by definition (5.6), implies (6.4) with u = A/8. □ 

7. Invariant measures. Let us consider the autonomous version of the 
system (5.1): 

dut = vtdt, 

(7.1) 

dvt = [Aut - 2avt + F{ut)] dt + S(nt) dWt, t > 0, 

with a given initial state {uq;vq), where F and S do not depend on t explic- 
itly. 

Let (pt = {ut;vt) and rewrite the system (7.1) as an evolution equation in 
the differential form, 

(7.2) d(l)t = A(l)tdt + J'{(bt)dt + dMt{(l)) 

with (j)Q = {uo;vo), where we set 



ut 
vt 



and 



A 





Fiu) 

■ 

A 



Mti 





Mtiu) 



I 

-2al 



where / is the identity operator on H. 

Under Conditions C given below, as an Ito equation in a Hilbert space, the 
solution (f)t, if it exists, is a Markov diffusion process in 7i (see [6], Chapter 
9). The transition probability function is given by 

Pt{^;B) = P{cl)teB\ct>o = 0, eeW, Bea{n). 

Suppose there exists an invariant measure fi on {7i,a{Ti.)), where cr(7i) de- 
notes the Borel fj-field of 7i. Then it satisfies ([7], page 12): 

fiiB)= [ Pti^;B)ii{dO yBGa{n). 
Jn 
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To show the existence of an invariant measure, we shaU speciahze Con- 
ditions A by assuming that the nonhnear terms satisfy a uniform Lipschitz 
continuity condition. To be precise, assume the fohowing Conditions C: 

(CI) Let -F(-) : — > H and S(-) : — > H, and let there exist positive con- 
stants bi,Ci for i = 1, 2, such that 

||i^('u)f <6i||n||?+ci 

and 

\\^t{u)\\l<b2\\u\\l + C2 

for any u E . 
(C2) There exist positive constants ki,k2 such that 

\\F{u) - F{u')\\^ < ki\\u-u'\\l 

and 

\\j:{u)-T,{u')\\'ji<k2\\u-u'\\l 

for any u, u' G H^. 
(C3) The constants 5j and ki satisfy 

{bi+b2X)A{ki+k2\)<Y- 

To show the existence of an invariant measure, we shall follow an approach 
by Da Prato and Zabczyk ([7], Theorem 6.3.2) for some stochastic dissipative 
systems. Though not directly applicable to the present problem, it can be 
adapted to proving the following theorem. 



Theorem 7.1 (Invariant measures). Suppose that the system (7.1) sat- 
isfies Conditions C. Then there exists a unique invariant measure fi on 
{TC,a{TC)). Moreover, given any bounded Lipschitz continuous function G 
on 7i, there are positive constants C and a2 such that 



(7.3) 

for any t> and ^ &7i. 



G{rj)Pt{i-dr,)- / G(r?)^(dr?) 
H JH 



<C(l + ||e||)e-°^* 



Proof. To extend the time domain for the system (7.1) to the whole 
real line R, introduce an independent copy Vt of the Wiener process Wt for 
t > 0. Define Wt by 



(7.4) 



W, 



V-u 



for t > 0, 
for t < 0, 
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and let Pt = cf{Ws : s < t} for t e R. Now, for t > r, let <l)t{T; = vt){T] = 
(nt(r; t't(T; ^)) be the solution of the extended system 

dut = vtdt, 

(7.5) dvt = [Aut-2avt + F{ut)]dt + J:{ut)dWt, t>T, 

where ^ = (Ci;^2) G 

Similarly to the derivation of the inequality (5.4) in the proof of Theorem 
5.1, it can be shown that, for A < {a A 7]i/Aa}, 



de^{ut;vt) < - Xe^ {ut; vt) dt + 



(7.6) 



^\\v^f + ^\\F{ut)f + mut)\\ 



dt 



+ 2{v^,j:iut)dWt), 



where dM^{u;v) = {v^,'E{ut)dWt). By making use of conditions (CI) and 
(C2), the above yields 



de^iut, vt) < —Xe^{ut; vt) dt + 
+ 2dMt{u;v) 



.A||2 



+ 



dt 



< 



dt + 2dM^{u;v), 



so that 



Ee\MsM < e^(6e-'(*-^)/' + (y + ^2) [ e-^^'^^^^ dr. 
Therefore, there exists a constant Ki > such that 



(7.7) Ee\(l)t{s,0]<Ki{l + e^{i)} for any t > s. 
For Ti>T2> 0, let 

u\=Ut{-Ti-C), vl = Vt{-Ti-C) for t>-T2, 

with i = l,2 and 

ut = uj-uf, dt = vl-v'^. 
Then it follows from (7.5) that we have 
dut =vtdt, 

(7.8) dvt = [Avt-2avt + 6F{ul;uf)]dt + 5T,{ul;uf)dWt, t> -T2, 
n-ra = {u\ - 6), ^-T2 = iv\ - 6), 
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where 
(7.9) 



Let v''^ = v + Xu and = e(nt; v^). As in (7.6), we can obtain the energy 
inequahty 



de^ < -Xe^dt + 



(7.10) 



^m'+\\muhum'+\mnhu',)\\i 



dt 



+ 2{i^,5J:{ut)dWt). 
In view of (7.9) and conditions (C2) and (C3), equation (7.10) yields 



de^ < -Xe^dt + 



X 



X. 



-,A||2 



+ 



\ut\\l 



dt + 2{v^,6^{ut)dWt) 



<--eUt + 2{v^,dJ:{ut)dWt), 



which imphes that 



In view of the bound (7.7) and the initial conditions in (7.8), we can show 
that 



Ee^ 



Ee\u\ - ulv] - vl) < 2K^{1 + e^(e)}e- 



A(i+r2)/2 



SO that 
(7.11) 
Let 

By setting t = in (7.11), we obtain, for any ^ G 

(7.12) Ee\^Pr, - V'n) < 2Ki{l + e^(0}e-^"^/', 

which goes to zero as r2 ^ oo. By Lemma 3.3 with A > 0, the energy func- 
tions e'^ and e are equivalent so that e(-) < Ce^(-) for some constant C > 0. 
Since e defines the energy norm || • || on by = e{(j)),(p G Ti, the set 
{tpT : r > 0} is a Cauchy family of random var iables in L^{n;n). Therefore, 
there exists a unique random variable ipoo G L'^{0,;TC) such that, for any 



lim E\\ 



0. 
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Since the random variables ipr = 0r(O;^) in distribution, ^,-(0;^) con- 
verges weakly to i/'oo as r ^ oo. It follows that Pt{^;-) converges weakly to 
the probability measure /x for ipoo as t — > oo, and fj, is the desired invariant 
measure for the system (7.1). 

To verify (7.3), let G : — > R be bounded and Lipschitz continuous on TC 
such that 

m)-G{v)\<lU-ri\\. 
Then, for any t > s > 0, we have 

/ G{ri)Pt{C;dr])- f G(7?)P,(e; dr?) 
Jn JH 

= \EG{i^t)-EG{^s)\ 

<K{l + e^(e)}^/'e-^^/^ 

due to the bound (7.12), for some constant K > 0. Therefore, by letting 
s — > t and invoking Lemma 3.3, inequality (7.3) follows. □ 

8. Examples. To illustrate the application of the stated theorems, let us 
specialize yl = (A — 1) in equation (4.1) to get 

d^u{x,t) = (A — l)u — 2adtu + f{u,x,t) + a{u, Du,x,t) dtW{x,t), 

0<t<T, xeVcRf^, d<3. 

(8.1) 

u{x, 0) = uo{x) eH^n L^", dtu{x, 0) = vo{x) G H, 
u{-,t)\dD = 0, 

where n > 1 . We shall present two examples with different kinds of nonlinear 
terms. 

Example 1. Let 

/(n, X, t) = -K-u^"-! + /3(x, t)u"', 

(8.2) 

a{u,Du,x,t) = C{x,t){l + \Du\'^)^u'', 

where k > is a constant, and (3 and C are some functions on D x as 
yet to be specified. The positive integers n, m and r are given such that 
2 < m < n, where n is any natural number for d <2 but, for d = 3, n <2. 
This is so because, according to a Sobolev lemma (see [3], Lemma 4.2), an 
L^-norm, depending on p and d, can be dominated by an iJ^-norm. Also we 
assume 5 G (0, ^) and < k < m(l — 26). 
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Owing to the Sobolev lemma mentioned above, as in [3], we can show 
that Conditions A are satisfied so that the equation has a unique global 
solution. In view of (8.2) and condition (A2), we set ^'{u) = 2Kn^"~^ so 
that ^{u) = ^IIm^IP and pt{u) = l3{-,t)u"^. Therefore, we have 

(8.3) (cI>'(n),u) = 2K||n"||2 and \\pt{u)f = m;t)u^f. 

By means of an elementary Young inequality ([9], page 61), it can be shown 
that, for any e > 0, we have 

,m\2 ^ ^„2n , /^^"^ 



qei/i' ' 

where q = n/{n — m) and q' = n/m, so that 

(8.4) \\Pt{u)f= I /3V'"dx<e||n"f + Ci / /^^^dx 

Jv Jv 

for some constant Ci > 0. Next consider the term 

||Sf(u)|||= / r{x,x)C\x,t){l + \Du\'^f^u^''dx. 
Jv 

By a repeated application of the Young inequality, we can deduce that, 
for any e',e" > 0, there exists C2 > such that 

(8.5) ||Sj(n)|||< / r{x,x)C^{x,t){e'u^'' + e"\Du\'^ + C2}dx. 

Jv 

Suppose that (3, C and r are bounded and continuous such that 

(8.6) \P{x,t)\<po,\C{x,t)\<Co and \r{x,x)\<ro 

for any x £T>,t>0. In view of (8.6), inequalities (8.4) and (8.5) yield 

\\pt{u)f<-^n) + Ci f P^i{x,t)dx, 
n Jv 

(8.7) 

||St(n)|||<roCo'|^^(^) + e"l|n||?| + C2ro^C'(a;,t)dx. 
From (8.3) and (8.4), using the notation in condition (B2), we see that 



f3i 


= 2k, 


71 


= 5i = 0, 


P2 


en 








72 


= 0, 




n 








2 e'k 
n 


73 





and 

(8.9) e{t) = Ci f P^'i{x,t)dx, p{t) = C2ro f C^{x,t)dx. 
Jv Jv 
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Therefore, condition (B3) takes the form 

2kA' - roCc 

i.W) 



2AtA2-roCo— A- — > „ 
n n 2 



roCoV'A<^. 

Since e, e' and e" are arbitrary, they can be chosen so small that condition 
(B3) holds simply for k > |. Assume this is the case. Then, by applying 
Theorems 4.1, 5.1 and 5.2, depending on the properties of 9 and p as defined 
in (8.9), we can draw the following conclusions: 

1. By the conditions in (8.6), it is clear that both 6 and p are bounded on R"*" 
so that, by invoking Theorem 5.1, we can conclude that the solution of 
the problem (8.1) is bounded in mean-square and there exists a constant 
Ki > such that 

snpE{\\ut\\l + \\dtUtf + \\u''f}<Ki. 

t>o 

2. Suppose that the functions 9 and p defined by (8.9) belong to -L^(R"'") so 
that 

/ / (3'^'^{x,t)dxdt<oo, / / (^{x,t)dxdt<oo. 
Jo Jv Jo Jv 

Then, by Theorem 5.2, the solution is ultimately bounded in mean-square 
if K > I and, furthermore, there is K2 > such that 

Esup{\\ut\\l + \\dtUtf + \\u''f}<K2. 

t>0 

3. Note that (8.1) has a null solution {u;v) = (0;0). Assume there exists 
ao > such that (e°o*6') and (e^^V) belong to L^(R+) or 

00 r roo r 

/ e°o*/32«(x,t)(ix(it <oo, / / e"o*C^(x,t)(ix(it < 00. 
Jv Jo Jv 

Then Theorem 6.1 shows that the null solution is exponentially stable in 
mean-square and, moreover, according to Theorem 6.3, the solution is in 
fact a.s. exponentially stable. 

Example 2. Consider a mildly nonlinear equation of the form 
d'iu{x, t) = {A- l)u - 2a dtu + /(n) 
^.11) +a{Du)dtW{x,t), 0<t<T, xeVcK^, 

u{-,t)\dD = 0, 
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subject to the initial conditions u{x, 0) = uo{x) G and dtu{x, 0) = vo{x) G 
H ^ where 

f{u) = — Kntan~^(l + u^), 

(8.12) 

(t{Du) dtW = ai{l + |i:)n|2}i/2 dtW^^^ + da dtW'^'^\ 

In the above equations k > 0, cii and (T2 are some constants, and W^^\-,t), 
are independent Wiener random fields with bounded, continuous 
covariant functions ?"i,r2, respectively. Rewriting (8.11) in the system form 
(7.1) and noting (8.12), it is easy to verify that 

(8.13) \\F{u)f<(^y\\u\\l 
and 



2viJ 

S(?x)|||j = erf / ri(x,x)(l + |L)iip) + fjf / r2{x,x)dx 
Jv Jv 



JV Jv 

for some C2 > 0, where r/i is the smallest eigenvalue of^ = (— A + 1) and 
ro = sup^jg-p ki(a;,a;)|. Similarly we can obtain the bounds 

\\F{u)-F{u')r<'^{l + l)'\\u-u'\\l 

(8.15) 

||S(Z)u) - S(Du')||| < ctItqWu - u'Wl 

for any u,u' G H^. In the notation of Conditions C, we can read off from 
(8.13) to (8.15) and find ci =0, 



2 



&1=(^) , ?^2 = CTiro, 



To satisfy condition (C3), we require that 



and 



2riiJ 2 



Then Theorem 7.1 (see the remark following Theorem 5.1) ensures the ex- 
istence of a unique invariant measure ^ in the state space TC for (8.11) and 
the the corresponding transition probability converges weakly to fi at an 
exponential rate. 
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APPENDIX 

Let W{x,t) be a continuous Wiener random field as given in Section 3. 
Then it may be regarded as an i?-valued Wiener process with a finite-trace 
covariance operator R with kernel r{x,y). We first define the stochastic 
integral with an a.s. bounded integrand. To this end, let a{x,t) be an a.s. 
bounded, continuous predictable random field such that 

(A.l) E / IWtfdt < oo. 

Jo 

We may consider at as a linear operator in H such that [crth](x) = (j{x, t)h{x) 
for any h€ H. Then it is easy to check that at ■ H ^ H is Hilbert-Schmidt 
a.s. and, noting (A.l), 

E [ Tr{atRa^) dt = E [ [ r{x,x)a^{x,t)dxdt 
Jo Jo Jv 

rT 

<roE / \\a^fdt< oo, 
Jo 

where * denotes the conjugation. Therefore, the stochastic integral 

(A.2) M{x,t)= f\{x,s)W{x,ds) 

Jo 

or 

Mt= f'a,dW, 



is well defined as a continuous //-valued martingale with mean zero and 
covariation operator Qt defined as (see [6], page 90) 

(A.3) {{M.,g),{M.,h))t= f\Qrg,h)dT, 

Jo 

where Qs has the kernel q{x,y, s) = r{x,y)a{x, s)a{y, s). 

Now we shall define a stochastic integral with an L^'-bounded integrand 
as shown in the proof of the following theorem. 

Theorem A.l. Let W{-,t) be a continuous Wiener random field with a 
bounded covariance function r{x,y) such that 

(A. 4) sup r(x, x) < To. 

Suppose that at = cr(-,t) is a predictable, continuous H -valued process sat- 
isfying the condition 

(A.S) E r \\a{;t)\\Pdt = E r [ \a{x,t)\Pdxdt<oo 

Jo Jo Jv 
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for an integer p> 2. Then the stochastic integral Mt in (A. 2) is well defined 
as a continuous H-valued, L'p -martingale with mean zero and covariation 
operator Qt for t S [0,T], as given by (A. 3). 

Proof. Since the set of bounded continuous functions on D is dense 
in L^iV) ([1], page 28), by smoothing, there exists a sequence {erf} of pre- 
dictable continuous random fields converging to at such that it satisfies 
condition (A.l) and 

(A.6) lim^E j'^ Wa'l -atfdt = Q. 

Therefore, as in (A. 2), the stochastic integral 

= f\^dWs 
Jo 

exists as a continuous if- valued martingale for each n. Let Mj. denote the 
Banach space of continuous L^-martingales Nt£H with norm ([6], page 79) 

i/p 

\N\\t = {E sup \\Nt\\P' 



L 0<t<T 

Then the sequence {M"} belongs to Aij-, since, by the B-D-G inequality, 
||M"||^ = ^ sup \\MJ}\\P 

0<t<T 

<CpE^^ J r{x,x)\a''{x,t)\'^dxdty^ 

<Cp{T)E r Wa'^fdt, 
Jo 

in which, in view of (A. 5), the upper limit is bounded, where Cp,Cp{T) are 
some positive constants. 
Now, for n> m, 

WW - M^'lfT = E sup ||Mj"-Mrir 

0<t<T 

T > p/2 



where 



TrQ^'^= / q'^''{x,x,s)dx 
Jd 
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JD 





It follows from (A. 4) and (A. 6) that 



M"" -M"'\\P.<Cp{T)E [ 




\Pdt 



Jo 



which goes to zero as n > m ^ oo due to (A. 6). Therefore, the sequence 
{M"} converges to the limit denoted by Mt, which is defined as a stochastic 
integral given by (A. 2). We can check that it preserves the properties of 
as stated in the theorem. □ 
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